Numerical study of velocity statistics in steady counterflow quantum turbulence 
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We investigate the velocity statistics by calculating the Biot-Savart velocity induced by vortex 
filaments in steady counterflow turbulence investigated in a previous study [Phys. Rev. B 81, 104511 
(2010)]. The probability density function (PDF) obeys a Gaussian distribution in the low- velocity 
region and a power-law distribution v~ 3 in the high-velocity region. This transition between the 
two distributions occur at the velocity characterized by the mean inter- vortex distance. Counterflow 
turbulence causes anisotropy of the vortex tangle, which leads to a difference in the PDF for the 
velocities perpendicular to and parallel to the counterflow. 

PACS numbers: 



Quantum turbulence, which is the disordered motion 
of a tangle of quantized vortices, was first observed in 
a counterflow experiment by Vinenii Since then, several 
properties of counterflow turbulence 2 - have been revealed 
through numerous experimental, theoretical, and numeri- 
cal studies. Although most counterflow turbulence stud- 
ies have focused on statistical values, such as the vor- 
tex line density L and the anisotropy of the vortex tan- 
gle, few studies have considered the PDF of the super- 
fluid velocity, which is important in classical turbulence. 
There are a number of reasons for this. The first reason 
is that the velocity field of quantum turbulence is not 
measurable experimentally. Since the superfiuid has no 
viscosity, particle image velocimetry (PIV), which is a 
powerful method for observing the velocity field, cannot 
accurately display the superfiuid velocity field. Another 
reason is that the numerical simulation of steady counter- 
flow turbulence has not yet been performed successfully. 
Schwarz performed the numerical simulation of counter- 
flow turbulence using a vortex filament model with lo- 
calized induction approximation, which neglects the in- 
teraction between vortices. As such, he was not able to 
obtain a statistically steady state in periodic boundary 
condition. 3 In a previous study 4 , we performed the nu- 
merical simulation of steady counterflow turbulence with 
the full Biot-Savart law, and the statistical values, such 
as vortex line density and the anisotropy, agree with the 
experimental observation. Thus, we can now numerically 
confirm the PDF of the superfiuid velocity field in steady 
counterflow turbulence. 

Recently, the velocity statistics of quantum turbulence 
has been discussed vigorously. Paoletti et aZ£ performed 
visualization of quantized vortices in a relaxation pro- 
cess of counterflow using solid hydrogen particles and ob- 
tained the non-classical (non-Gaussian) PDF of the par- 
ticle velocity. They reported that the non-classical statis- 
tics are due to the velocity induced by the reconnection 
of a quantized vortex because the PDF exhibits a power- 
law distribution of v~ 3 , which is derived from the vortex 
velocity before or after reconnection. However, they ob- 
served the velocity of particles, which is not necessary 
the velocity of the superflow. The non-classical velocity 
statistics were also confirmed by White et al& They per- 



formed numerical simulations of quantum turbulence in 
a trapped Bose-Einstein condensate by calculating the 
Gross-Pitaevski equation to obtain the PDFs of the su- 
perflow field. The PDFs show not classical Gaussian dis- 
tributions, but rather power-law distributions, which is 
due to the velocity field v = «/ (2nr) induced by the sin- 
gular quantized vortex, where k is the quantum of circu- 
lation and r is the distance from the core of a quantized 
vortex. 

In the steady state of counterflow turbulence, we show 
that the PDF obeys non-classical distributions due to the 
singularity of the quantized vortex. We also reveal that 
the PDF depends on a line density and the anisotropy of 
a vortex tangle. 

Counterflow is the internal convection between the su- 
perfiuid and the normal fluid having a relative velocity 
of v ns = v n — v sa , where v„ and v sa are the normal 
and superfiuid velocities generated by the counterflow. 
Hence, in the steady counterflow turbulence, the super- 
flow velocity produced at a point r by a filament is given 
by 
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where Si refers to a point on the filament and the inte- 
gration is performed along the filament. The first term 
denotes the Biot-Savart velocity induced by the fila- 
ment. We can obtain the superfiuid velocity field by 
calculating Eq. (1) from the configuration of the vor- 
tices. In the present paper, the steady counterflow tur- 
bulence is calculated in the same manner as in Refer- 
ence 4, in which we performed numerical simulations of 
steady counterflow turbulence using the full Biot-Savart 
vortex filament model, and the statistical values (vor- 
tex line density and anisotropy) agreed quantitatively 
with the experimental observations. In order to inves- 
tigate how the PDF depends on the vortex line den- 
sity, we calculate the velocity field under two conditions, 
namely, v sa = (0, 0, -0.496)cm/s (L ~ 10, 000cm~ 2 ) and 
v sa = (0,0, -0.310)cm/s (L ~ 4200cm- 2 ) at T = 2. IK. 
The vortex configuration at the steady state is shown in 
FiglU The vortex tangles are quite anisotropic in the 
direction of v ns = v„ — v sa (z direction). Therefore, we 
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also investigate the influence of anisotropy on the PDF. 




FIG. 1: (Color online). Configurations of the steady state of 
counterflow turbulence at T = 2. IK. The counterflow relative 
velocity w sa is in the direction of —z. (a) v 3a = —0.310 cm/s, 
and (b) v sa = —0.496 cm/s. 



We calculate the superflow velocity field in steady 
counterflow turbulence by the spatial discretization on 
a 100 3 grid. Figure [2] shows the PDFs of v x , v y , and v z . 
Although it is known that the PDF of classical viscous 
turbulence is Gaussian^^ the PDFs of the present coun- 
terflow calculation exhibit a non-Gaussian distribution 
with a large tail in the high- velocity region. Since the vor- 
tex tangle of steady counterflow turbulence is isotropic in 
the direction perpendicular to the relative velocity v ns , 
the PDF of v x almost overlaps with that of v y , with the 
peaks of two PDFs at v x = and v y — 0. In contrast, 
since the superfluid velocity v sa — —0.496 cm/s due to 
counterflow is applied in the —z direction, the PDF of v z 
has a peak at v sa . 
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FIG. 2: (Color online). Probability distribution of a velocity 
component v x , v y , and v z in steady counterflow turbulence of 
v sa = —0.496 cm/s at T = 2. IK. The vertical dot-dashed line 
indicates v sa - 

Figure [3] shows the PDFs inverting the negative re- 
gion around v x — v v = and v z o = v z — v sa = cm/s. 
Since the vortex tangle is isotropic in the direction per- 
pendicular to v ns , the PDF of v x is symmetric to v x = 0. 
The PDF of v z is asymmetric to the peak v z q — be- 
cause the counterflow turbulence is anisotropic in v ns 
direction. In this anisotropic state, many vortices lie in 
the plane perpendicular to v ns . In counterflow, mutual 



friction can expand only the vortices for which the self- 
induced velocity is opposite in direction to v ns , so that a 
greater number of such vortices survive, as compared to 
vortices for which the self-induced velocity is parallel to 
v ns . Hence, the velocity field of the counterflow vortex 
tangle is anisotropic in the z direction. 
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FIG. 3: (Color online). Probability distributions of (a) v x 
and (b) v z , (v sa — —0.496 cm/s, T = 2.1 K) inverting the 
negative region around (a) v x = and (b) v z o = cm/s. 

Figure [4] shows double logarithmic plotted PDFs of 
the v x and v z components. Both distributions exhibit a 
power-law distribution of Pr(^i) oc v^ 3 (i = x, y, zO) in 
the higher- velocity region. These velocity statistics have 
been investigated in previous studies.— For the single, 
straight vortex case, the probability of separation occur- 
ring between r and r + dr is 2irrdr, and the velocity scales 
as 1/r, which leads to Pr(v) ~ Pr(r(v))\dr / dv\ ~ 1/v 3 . 
The PDF converges to a Gaussian distribution in the 
low-velocity region, probably because the vortex config- 
uration is random in the tangle. This property of the 
PDF was investigated in a previous studyi^ in which the 
authors numerically set N vortex points at random posi- 
tions, and the PDF of the velocity at the vortex point was 
obtained. The transition of the PDF between a Gaussian 
distribution and a power-law distribution of v~ 3 was also 
confirmed. 

We can roughly estimate the transition velocity from 
the Gaussian distribution to the power-law distribution. 
In order to easily understand the velocity field induced 
by multiple vortices, we consider the simple case of two 
straight parallel vortices, as shown in Fig. [5] Although 
primarily the 1/r velocity appears near each vortex, 
in the halfway range between vortices, the velocity be- 
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comes complicated because the velocities induced by the 
two vortices become comparable and interfere with each 
other. Hence, the statistics of velocity appear to change 
near the midpoint between vortices. In the vortex tangle, 
the mean inter- vortex distance is denoted by I — 1/L 1 / 2 , 
and so the midpoint between the vortices is located at 
1/2. Thus, the transition velocity of the statistics should 
be represented by 



In Fig. HI the vertical solid lines refer to the estimated 
transition velocity Vt- The actual transition velocity for 
the power-law distribution approximately agrees with the 
estimated transition velocity v t - The mean inter- vortex 
distance in the z direction is different from that in the 
x and y directions due to the anisotropic vortex tangle, 
which leads to a difference in the actual transition veloc- 



ity between the PDFs for v x and v z . 

In conclusion, we investigate the velocity statistics by 
calculating the Biot-Savart velocity induced by vortex 
filaments in steady counterflow turbulence investigated 
in a previous study^ The probability density function 
(PDF) obeys a Gaussian distribution in the low- velocity 
region and obeys a power-law distribution u -3 in the 
high-velocity region. The transition between these two 
distributions occurs at the velocity characterized by the 
mean inter- vortex distance. The anisotropy of the vor- 
tex tangle due to counterflow leads to a difference in the 
PDF between the velocities perpendicular to and paral- 
lel to the counterflow. In the future, we intend to study 
how the velocity statistics are related to other statistical 
quantities, such as energy spectral 
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FIG. 4: (Color online). Log- log plot of the probability distri- 
butions of (a) v x and (b) v z q in the steady state vortex tangle 
calculated under v sa = —0.496 cm/s (L — 10, 000 cm -2 ) and 
v sa = -0.310 cm/s (L ~ 4200 cm" 2 ). The blue and red 
solid lines indicate the vt derived from the vortex line den- 
sity L = 10000 cm -2 and L = 4200 cm" 2 , respectively. The 
dashed lines show oc v~ 3 , and the dot-dashed lines show the 
Gaussian distribution. 
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FIG. 5: (Color online). Schematic diagram of the velocity 
field induced by two parallel vortices. The color density cor- 
responds to the magnitude of the velocity induced by the 
vortex, which attenuates with 1/r. 



